Mon. Not. R. Astron. Soc. 000, [[HS] (2012) Printed 1 January 2013 (MN JATfcjX style file v2.2) 



Rotational threshold in global numerical dynamo 
simulations 



(N 

o 

(N 

O 
0) 

Q 



m 

6 

o3 



> 

O 

On 

(N 



13 



M. Schrinner 1 * 

1 MAG(ENS/IPGP),LRA, Ecole Normale Superieure, 24 Rue Lhomcmd, 75252 Paris Cedex 05, France 



Accepted . Received; in original form 



ABSTRACT 

Magnetic field observations of low-mass stars reveal an increase of magnetic activ- 
ity with increasing rotation rate. The so-called activity-rotation relation is usually 
attributed to changes in the underlying dynamo processes generating the magnetic 
field. We examine the dependence of the field strength on rotation in global numerical 
dynamo models and interpret our results o n the basis of energy c o nside rations. In 
agreement with the scaling law proposed by IChristensen fc Auber t (2006), the field 
strength in our simulations is set by the fraction of the available power used for the 
magnetic field generation. This is controlled by the dynamo efficiency calculated as 
the ratio of Ohmic to total dissipation in our models. The dynamo efficiency grows 
strongly with increasing rotation rate at a Rossby number of 0.1 until it reaches its 
upper bound of one and becomes independent of rotation. This gain in efficiency is 
related to the strong rotational dependence of the mean electromotive force in this 
parameter regime. For multipolar models at Rossby numbers clearly larger than 0.1, 
on the other hand, we do not find a systematic dependence of the field strength on 
rotation. Whether the enhancement of the dynamo efficiency found in our dipolar 
models explains the observed activity-rotation relation needs to be further assessed. 
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1 INTRODUCTION 

There is considerable observational evidence that magnetic 
activity on low mass stars increases with increasing stel- 
lar rotation rate until it saturates an d reaches a co nstant 
level for very fast rotating stars (e.g. iReinersI I2012I ). This 
so-called activity-rotation relation is primarily based on 
observations of chromospheric or coronal magnetic activ- 
ity indicators, i.e . on chromospheric or coronal emission 
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i Skumanichl Il972l; iNoves et all Il98l ; iDelfosse et al. 1 1 19981 ; 
iPizzolato et al.|[2003l ). Furthermore, it is supported by mag- 
netic flux measurements in stars of spect r al typ es G-M (|Saarl 
l200ll ; iReiners et aHl2009h . INoves et all (| 19841 ) pointed out 
that chromospheric emission and thus magnetic activity de- 
pending on spectral type and rotation is well described by 
a single parameter, the Rossby number Ro = P ro t/T c , with 
P ro t being the observed rotation period of a star and r c a 
convective turnover time derived from mixing-length theory. 
The observed increase of magnetic activity with decreas- 
ing Rossby number is usually attributed to changes in the 
underlying dynamo processes generating t he magnetic field 
i|Donati fc LandstreetJ 120091 ; IReiners! |2012| ). However, until 
now, theoretical arguments explaining the activity-rotation 
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relation are poorly developed and formulated only on a 
heuristic level. 

The standard argument worked out in the framework of 
linear mean-field theory refers to the so-called dynamo num- 
ber, D, a dimensionless parameter which determines the dis- 
tance to the dynamo threshold of a given dynamo model. If 
D exceeds some critical value, D > D c , a small initial mag- 
netic perturbation may grow exponentially, i.e. the field-free 
state is linearly unstable. The larger D, the larger are the 
growth rates of the magnetic field expected in this scenario. 
Under some simplifying assumptions, D scales i nversely pro- 
porti onal to the square of the Rossby number (|Noves et al.l 
Il984h . Therefore, one might argue that dynamo action can 
be easier excited and thus leads to larger field strengths if 
the Rossby number decreases. 

The principal objection against this argument is that it 
is entirely based on linear, kinematic theory. For stellar dy- 
namos with D ^> D c , the magnetic field would grow rapidly 
until the Lorentz force changes the velocity field and as a re- 
sult the magnetic field saturates. In this dynamical regime, 
linear theory is in general no longer applicable. Even if the 
dynamo number predicted the onset of dynamo action prop- 
erly, its predictive power for the saturation level of the mag- 
netic field would remain uncertain. 

A successful scaling law for the field strength of fast 
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rot ating stars and pla n ets, o n the other hand, was presented 
by IChristensen et aTl (j2009t ) . They showed that a scaling 
originally deriv ed from geodynamo models in the Boussinesq 
approximation (IChristensen fc Aubert|[2006l ) also applies to 
certain classes of stars. The scaling law reads 



(1) 



and is based on the available energy flux q c , which sets 
the field strength apart from the density, g, and a convec- 
tive length scale relative to the temperature scale height, 
L/Ht- The coefficient / h m ^ f gives the ratio of Ohmic 
to total dissipation and was set to unity in this context. 
IChristensen et aD (|2009l ) highlighted two remarkable find- 
ings: First, the flux based scaling law is at first glance inde- 
pendent of rotation, and second, it is only applicable to fast 
rotating stars and planets. 

In the following we will argue that the flux based scal- 
ing law (0 is in fact valid for a wide range of rotation rates, 
but it is not independent of rotation. On the contrary, we 
claim that the rotational dependence is captured by / h m 
and simply eliminated by setting / h m to unity. This in- 
deed seems to be justified only for stars in the rotationally 
saturated regime. Our analysis is based on 30 numerical dy- 
namo models in the Boussinesq approximation. The model- 
ing strategy and the models are briefly described in the next 
section. Hereafter, we present results for / h m revealing its 
dependence on rotation rate and discuss the implications of 
our finding for the activity-rotation relation. 



2 DYNAMO CALCULATIONS 

Our dynamo models are solutions of the MHD-equations for 
a conducting Boussinesq fluid in a rotating spherical shell. 
Convection is driven by an applied temperature difference 
AT between the inner boundary at radius Ti and the outer 
boundary at r . The governing equations for the velocity v, 
the magnetic field B, and the t emperature T writ ten in a 
dimensionless form proposed bv lOlson et al.l (|l999h are 

E (J^- + v ■ Vd - V 2 u^ + 2z x v + VP = 

Ra— T+ -!-(V > 

r Pm 

dT 

V x (v x B) + 
V • v = 0, 



dB 
~dt 



< B)x B , 


(2) 
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(3) 




(4) 


V • B = 0. 


(5) 



They are controlled by four dimensionless parameters, the 
Ekman number E = v/flL 2 , the (modified) Rayleigh num- 
ber Ra = argoATL/vQ, the Prandtl number Pr — v/k, 
and the magnetic Prandtl number Pm — f/rj. In these defi- 
nitions, L denotes the width of the convection zone, Q stands 
for the angular velocity, a is the thermal expansion coeffi- 
cient, g is the gravitational acceleration at the outer bound- 
ary, and v, rj, k are the kinematic viscosity, the magnetic and 
the thermal diffusivity. The mechanical boundary conditions 
are no slip and the magnetic field continues as a potential 
field outside the fluid shell. 

In the following paragraph we need to define some fur- 
ther quantities which are used throughout in the paper. The 



time-averaged ratio of ohmic to total dissipation, /ohm, is 
computed as 

Wj_ 

W b 



John 
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with the rate of ohmic dissipation 

^ = ^L_/(VxS) 2 dv, (7) 
and the power Wb generated by buoyancy forces, 

Wb = / -VrTdv. (8) 

E J T 

Moreover, we define a Rossby number for our models 
similar to the observational one. It is given by the ratio of 
the Rossby radius to a typical convective length scale l c , 
Roe = D rms / '(Q£ c ), where v ims stands for the rms velocity 
of the flow and £ c is derived from the kinetic energy sp ec- 
trum (|Schrinner et al.ll2012l ; IChristensen fc Aubertteooa ). A 
non-dimensional measure for the convective energy flux is 
the Nusselt number, Nu, defined as the ratio of the total 
heat flow to the conducted heat flow. Finally, we note that 
the temperature scale height Ht for our models is given by 
Ht = cp/(ag ) with the heat capacity cp. 

We aim at comparing models which differ only in their 
rotation rates. Therefore, we keep the thermodynamically 
available energy flux and the diffusivities for a sequence of 
models constant and vary only the angular velocity. Trans- 
lated in non-dimensional quantities, this means: We keep (in 
a first approximation) the Rayleigh number over some criti- 
cal Rayleigh number Ra c , the Prandtl number and the mag- 
netic Prandtl number constant and change successively the 
Ekman number. More precisely, we try to keep the Nusselt 
number Nu for a sequence of models constant. Because Nu 
is an output parameter in our simulations and only roughly 
determined by Ra normalized by its critical value, Ra/Ra c 
has to be adjusted accordingly. Some more details are given 
in the Appendix. 

Due to computational limitations all current numerical 
dynamo simulations run in a parameter regime which is not 
appropriate for stellar interiors. Moreover, Boussinesq mod- 
els do not account for the strong density variation in stars 
and thus certainly do not reproduce stellar dynamo pro- 
cesses in realistic detail. However, our models are adequate 
to study the flux-based scaling law which was originally de- 
rived from Boussinesq models, too. 



3 RESULTS 

We considered sequences of models with Nu ~ 2.2, 3.5, and 
7 and Pm varying between 3 and 7. The Prandtl number 
Pr was always set to unity in our simulations. For given 
Pm and Nu, we obtained a sequence of models by varying 
the Ekman number between E — ICP 3 and E — I0 -5 ; some 
models with E = 3 • I0" 3 and E = 3 ■ 1CT 6 could also be 
added to our sample. On the other hand, simulations with 
Nu w 7 and E ^ 3 x 10 were numerically not feasible. The 
magnetic Reynolds number of our models, Rm = « rms L/r/, 
is always larger than 100 and thus far above the minimum 
value of Rm = 40 needed to obtai n dynamo action in this 
setting (|Qlson k, Christe nscn 2006j). We show in Appendix 
A that there is a slight but systematic increase of Rm with 
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Figure 1. ^°/( c / omI1 ) versus the flux based Rayleigh number. 

The constant c was chosen to be 0.92 for dipo l ar and 0.48 for mu l- Figure 2. The ratio of ohmic to total dissipation, / ohm , versus 

tipolar models l lChristensen fc Aube rq |2006| ; I Christensen! |2010|) . the Rossby number for sequences of models with Nu = 2.2 and 

The solid line is not a fit but the prediction of equation d 1 D . various Pm. 



rotation rate for a sequence of models at constant Nusselt 

and Prandtl nu mbers. 

Following IChristensen fc Auberd |2006) a non- 
dimensional form of the flux-based scaling law is obtained 
by dividing relation Q by gQ 2 L 2 , 

B 2 /(2fi gQ 2 L 2 ) ~ / ohm ( 9c /(^ 3 L 3 ) L/H T ) {2/3) . (9) 

With the Lorentz number Lo = B /(figQ, 2 L 2 ) 1//2 , the flux 
based Rayleigh number Rciq = q c ag /(4ivQCpQ. 3 L 2 ), and 
L/Ht = Lag /cp, relation Q may simply be written as 

Lo/fUl = cRa^ 3 . (10) 

IChristensen fc Auberd (2006) found the prefactor c to be 
0.92 for models with a predominantly dipolar field geom- 
etry and a s lightl y lower value of c = 0.48 was given by 
IChristensen! l|2010t ) for multipolar models. In Fig.[T]we plot- 
ted Lo/^cf^l^) against Rciq in logarithmic scales for our 
sample of models. Independent of their Ekman number, all 
models are in agreement with the scaling (|10|) indicated in 
Fig. [1] by the solid line. Apparently, the flux-based scaling 
law holds for all models independently of their rotation rates. 

However, the ratio of ohmic to total dissipation, /ohm, 
may increase drastically with rotation rate, as demonstrated 
in Fig. [2] Shown is / h m versus the Rossby number for se- 
quences of models with Nu = 2.2 and various Pm. At 
Rot ~ 0.12, the rate of ohmic diffusion increases rapidly 
until the steep slope flattens and /ohm saturates for lower 
Rossby numbers. On the other hand, a strong dependence 
of /ohm on Pm cannot be inferred from Fig. [2] 

The strong increase of / h m at Roi < 0.12 coincides 
with a transition from multipolar dynamo models at higher 
Rossby number to models with a dipole dominated mag- 
netic field. In fact, the dependence of /„i lm on the rotation 
rate changes crucially at the regime boundary as shown 
in Fig. [3] For a sequence of predominantly dipolar mod- 
els with a slightly higher Nusselt number, Nu — 3.5, the 
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Figure 3. The dynamo efficiency /ohm versus the Rossby number 
for sequences of models with Pm = 3 and Nusselt numbers Nu = 
2.2, Nu = 3.5, and Nu = 7. The vertical dashed line indicates 
the boundary between the dipolar (Roi < 0.12) and multipolar 
dynamo regime. 



increase of / h m is qualitatively reproduced. However, the 
amplified convective energy flux shifts the sequence towards 
higher Rossby numbers. For an even higher Nusselt number, 
Nu — 7, the sequence of models falls entirely in the multi- 
polar regime and a systematic increase of / D hm with rotation 
rate is no longer observed. 



4 M. Schrinner 



4 DISCUSSION AND CONCLUSIONS 

In an equilibrium state, the energy released by buoyancy in 
our models is dissipated by viscous dissipation and ohmic 
diffusion. The latter requires that a magnetic field is built 
up by dynamo action and the rate of ohmic dissipation de- 
termines the fraction of the available power used for the 
magnetic field generation. For fast rotators fohm increases, 
this means that a larger fraction of the available power is 
converted to magnetic energy and dynamo action becomes 
more efficient. According to relation fl]l. the growth of /ohm 
visible in Fig.[2]leads to an increase of the average magnetic 
field strength by an order of magnitude. Because /ohm is 
bound by one, the field strength saturates for even higher 
rotation rates and then becomes independent of rotation. 
Both, the increase at Roi « 0.1 and the saturation of the 
field strength are in good agreement with observations. Yet, 
not much is known about the saturation level of the mag- 
netic field in slowly rotating stars, except that it falls below 
the v alue predicted by flj with /ohm = 1 (fchr stcnscn et al 
2009). With a Rossby number of Rot > 0.5 (jReiners et al 
2012J), the Sun is an example for a slow rotator 
ing an energy flux of q c = 63MW/m 2 
p = 1.4kg/m 3 and an av erage internal field of B = 0.063 T 
l|Christensen et al1l2009l ). relation JT} requires / h m ~ 0.07. 
This would be consistent with the range of /ohm presented 
in this study. We note, however, that the flux based scaling 
law is som ewhat at odds with the estimate of fohm = 10 -3 
derived bv lRem pcl (2006) from dynamic flux-transport solar 
dynamo models. 

The decline of /ohm at Rot « 0.1 visible in Fig. [2] and 
Fig. [3]is related to a rotational dynamo threshold. It is char- 
acterized by a minimum magnetic Prandtl number, Pm cr it, 
below which self-sustain ed dynamo action does not occur. 
IChristensen et all l|l999h found that Pm cr it is a function of 
only the Ekman number and varies in the dipolar regime as 



Assum- 
a mean density of 



Pm clit = 450 E u 



(11) 



Models with given diffusivities approach this dynamo 
threshold if their rotation rate is decreased and /ohm drops 
to zero. 

Equation (|11|) is independent of any velocity ampli- 
tude and does not relate the rotational threshold to a given 
Rossby number. Hence, particular low values for /ohm could 
in principle be found at any Rossby number (and at any Rm 
). However, we could not confirm equation (|11|) for the mul- 
tipolar dynamo regime where /ohm remains low and does not 
change significantly with rotation rate (see Fig. [3]). Attempts 
to identify a similar decline of /ohm close to a rotational 
threshold at much lower Rossby numbers also failed. Mod- 
els in this parameter regime with a lower Ekman and mag- 
netic Prandtl often bifurcate subcritically ()Morin fc Dormvl 
2009). Consequently, the magnetic field strength and /ohm 
remain high and collapse abruptly to zero only beyond the 
dynamo threshold. 

What explains the rotational threshold and the high 
rotational sensitivity of fohm for models in the dipolar 
regime close to Roe ~ 0.1? Dipolar models typically ex- 
hibit different characteristic length scales for their veloc- 
ity and their magnetic field with £v < £b- The balance of 
the advection and the diffusion term in the induction equa- 
tion then leads to a modified magnetic Reynolds number, 



Rm = v rlnB £% / (r/ £y) , which needs to exceed a critical value 
for the onset of dynamo actio n. For given diffusivities, £% is 
inversely proportional to n rms ()Christensen fc Tilgnerll2004T ) . 
Therefore, an increase of the velocity amplitude is compen- 
sated by smaller £b and does not change Rm. This heuristic 
argument might explain why equation (lll[) is independent of 
any velocity amplitude and holds only in the dipolar regime. 
The strong dependence of the dynamo efficiency on rotation 
rate, however, requires some further explanation. 

Dynamo models in the dipolar regime with Roi < 0.1 
may be ade quately described in the framework of mean- 
field theory (IKrause fc Radlerl [TQ80I: ISchrinner et al.l [20071 ; 
ISchrinnerll201ll ; ISchrinner et al.ll201lr ). The mean-field for- 
malism provides useful concepts to better understand the 
influence of rotation on the dynamo processes in our models. 
It is usually set up by splitting the velocity and the magnetic 
field in a mean and a residual component varying on different 
length scales, v = V + v' and B — B + b' . Mean quantities 
denoted here by an overbar may be thought of as combined 
azimuthal and time averages <|Schrinnerll201ll) . They vary on 
a scale similar to the system size L. Fluctuating quantities, 
on the other hand, vary on the scale of the velocity field, £v- 
We emphasize that this is not an assumption but a result 
obtained from direct numerical simulations in a particular 
parameter regime. The mean flow V is negligible for the 
mode ls considered here ( (Olson et al.l 1 19991 ; ISchrinner et al.l 
120071 ). Therefore, the induction equation ((2| separated for 
the mean and the residual component may be written as 

d~B 1 

~~dt ~ ~Pm 

db' „ , L ,w 1 
A — V x (v x b ) - — — 

dt K ' Pm 



V 2 B = V x v x b' (12) 
' = V x (v x B). (13) 



The residual magnetic field b' is diffusive on the length scale 
£v and would rapidly decay without the source term on the 
right-hand side of l)13p. Thus, the magnetic field-generation 
depends decisively on B and in particular on the so-called 
mean electromotive force, £ = v x b' . Similarly, the mag- 
netic energy density is dominated by the mean field. The 
energy equation for B reads 



d_ 

dt , 



J 00 2 J y 



Edv, 



(14) 



where E is the mean electrical field, j is the mean current 
density and V denotes the volume of the fluid shell. With 
Ohm's law, 



j = Pm(E + £) 



equation (|14[) yields 



1 

~Pm 



j 2 dv = / j ■ £ dv. 
v Jv 



(15) 



(16) 



for an equilibrium state. Hence, also the mean Ohmic diffu- 
sion is controlled by the electromotive force, i.e. by the cor- 
relation of the residual velocity and the residual magnetic 
field, and thus linked to rotation. It is expected that rota- 
tion strengthens the correlation between v and b' . Indeed, 
£ grows with rotation rate for a sequence of models with 
fixed Nu, as shown in Fig. [4] Therefore, also Wj and even- 
tually the dynamo efficiency /ohm increase with decreasing 
Rossby number. We note, however, that also Wb (in units of 
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Figure 4. Left: Time-averaged rms values of the electromotive 
force in units of [v/L (griuQ) 1 / 2 ] for a sequence of models with 
Nu 2.2 and Pm = 4 versus the Rossby number. Right: The 
rate of Ohmic dissipation (stars) and total dissipation (squares) 
in units of gu 3 /L versus Roi for the same sequence of models. 



gv 3 /L) increases for this sequence of models, though some- 
what slower than Wj. 

For clarification, we stress that the decrease of the dy- 
namo efficiency with increasing Rossby number in our mod- 
els is not caused by an J?m-dependent quenching of the 
electromotive force, which is sometimes called catastrophic 

I n — II — I 

quenching (see Brandenb urg fc Su bramanian 2005, and ref- 
erences therein). In contrast to the catastrophic quenching 
scenario the mean electromotive force increases with Rm in 
our simulations (see also Appendix A). 

In summary, the field strength of our models is set by 
the available energy flux and via / h m by the rotation rate. 
The dynamo efficiency / h m increases strongly with rota- 
tion rate at Roe ~ 0.1 and saturates at smaller Rossby 
numbers. The high rotational sensitivity of / h m is related 
to a rotational dynamo threshold and finally to the strong 
dependence of the mean electromotive force on rotation in 
this parameter regime. For multipolar dynamos at higher 
Rossby number, however, the dynamo efficiency seems to 
be almost independent of rotation. Similarities with the ob- 
served activity rotation relation are encouraging and need 
to be further assessed. 
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APPENDIX A: SCALING OF Ra/Ra c AND Rm 
AT CONSTANT NUSSELT NUMBER 

We use scaling laws given by IChristensen fc Aubertl (|2006t ) 
to show that the ratio Ra/Ra c decreases slightly with rota- 
tion rate for a sequence of models at constant Nusselt and 
Prandtl numbers, whereas the magnetic Reynolds number 
increases. 

The Nusselt numb er scaling proposed by 
IChristensen fc Aubertl (|2006t ) may be written as 

Nu-1 ~ RaE, (Al) 

provided that Nu > 1 and convection is sufficiently super- 
critical. Moreo ver, the cri tical Rayleigh number varies as 
Ra c ~ E~ 4/3 (|Busselll97dh , and we finally find Ra/Ra c ~ 
E 1 ^ 3 for models at constant Nusselt number. For Nu ~ 2.2, 
3.5, and 7, we considered a maximum Rayleigh number of 

6, 15, and 50 times its critical value. 

The Rossby number scaling from lChristensen fc Aubertl 
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(2006) together with the Nusselt number scaling yields 

Ro ~ ((Nu ~ 1) E/Prf 77 . (A2) 

With Rm = RoPm/E we conclude that the magnetic 
Reynolds number increases slightly with increasing rotation 
rate, Rm ~ E' 023 . For iVu w 2.2 and Pm = 4, for in- 
stance, we obtained -Rm fa 130 at E = 10 -3 and 7?m « 270 
at £ = 10" 5 . 



